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ABSTRACT 



Sensitivity limits of ground-based infrared interferometers using aperture synthesis 
are presented. The motivation of this analysis is to compare an interferometer 
composed of multiple large telescopes and a single giant telescope with adaptive optics. 
K*- ' In deriving these limits, perfect wavefront correction by adaptive optics and perfect 

cophasing by fringe tracking are assumed. We consider the case in which n beams are 



CN ■ pairwise combined at n{n — l)/2 detectors ("C2 interferometer) and the case in which 



all the n beams are combined at a single detector ("C^ interferometer). Our analysis 
covers broadband observations by considering spectral dispersion of interference 
fringes. In the read noise limit, the "C„ interferometer with one dimensional baseline 
configuration is superior to the "C2 interferometer, while in the background limit, 
the advantage of the one-dimensional "C„ interferometer is small. As a case study, 
^H ' we compare the point-source sensitivities of interferometers composed of nine 10-m 

^ • diameter telescopes and a 30-m diameter single telescope with adaptive optics between 

1 and 10 /im for lO-o" detection in one hour. At J and H, the sensitivities of the 



interferometers are limited to 25 to 24 mag by read noise and OH airglow background, 
while that of the single telescope is limited to 28 to 26 mag by OH airglow background. 
C^ I Longward of 2 fiui, the sensitivities of the interferometers and the single telescope 

are all limited by instrumental thermal background. At K, the sensitivities of the 
interferometers are around 23 mag, while that of the single telescope is 27 mag. At N, 
the sensitivities of the interferometers are around 10.7 mag, while that of the single 
telescope is 14.4 mag. 

Subject headings: techniques: interferometric 



1. Introduction 

Fundamental limitations of optical and infrared aperture synthesis imaging had been an 
academic subject in which astronomers were not keenly interested. The analyses of sensitivities 



are mathematically rather involved and they appeared in journals of optics ( prasad fc Kulkarni 



(1989)1 (PK89), [Nakajima fc Matsuhara (2001)| (NMOl)) instead of those in astronomy. The 



analyses were academic partly because they dealt with ideal conditions in which atmospheric 
disturbances did not exist. Therefore their results were applicable only to space interferometers 
which were instruments in the distant future. They were also academic partly because optical and 
infrared interferometers were not considered as routinely used astronomical instruments unlike 
radio interferometers in radio astronomy. 



However, as pointed out by Roddier &: Ridgway (1999) (hereafter RR99), the situation has 



changed due to the completion of large ground-based telescopes with 8-10 m diameters and to the 
planning of next-generation giant telescopes with 30-100 m diameters. In the plans of the giant 
telescopes, they are all equipped with adaptive optics which deliver diffraction limited images. The 
giant telescopes are always useful as optical buckets. However for the purpose of high resolution 
imaging, it is not clear whether the synthesized images from an interferometric array of large 
telescopes or diffraction limited images of giant telescopes are superior. In terms of resolution, 
the interferometric array is obviously advantageous, but in terms of sensitivity or dynamic range 
the quantitative comparison is more complex. Optical/infrared interferometric arrays with small 
telescopes currently in operation or being built are dedicated to stellar astronomy. With the 
advent of an interferometric array composed of large telescopes, one essential question is whether 
interferometry finally has an impact on extragalactic astronomy. This question is directly related 
to the limiting magnitude. 

The development of the fringe-tracking technique has enabled to cophase a number of telescopes 
and a ground-based interferometric array can be regarded as an ideal interferometer without 
atmospheric disturbances. Therefore the sensitivity analyses made for space interferometers in the 
past can now be applied to ground-based interferometers as well. Under these circumstances, it is 
worthwhile to highlight the results of previous analyses of optical and infrared synthesis imaging 
in an astronomical journal and apply them to ground-based interferometry composed of large 
telescopes. 

RR99 pioneered the analysis of ground-based interferometry from the point of view of 
comparison with future giant telescopes. However their analysis is not of aperture synthesis and 
as we will show later that their results are not necessarily in agreement with ours. Although this 
paper mostly utilizes previously obtained formulae of signal-to-noise ratios (SNRs) in PK89 and 
NMOl, we also present new results to cover the full parameter space. 

The beam combination geometry is a major issue in studying the sensitivity of an 
interferometer. One extreme is the "C2 interferometer, in which n beams are divided into n{n — 1) 
subbeams and they are pairwise combined at "C2 = n{n — l)/2 detectors. There is one detector for 
each baseline. The other extreme is the "Cn interferometer in which all the n beams are combined 
at "'C„ = 1 detector. In this paper, we present the sensitivities of both interferometers. We 
consider idealized interferometers, but the beam-combination scenario for actual interferometers 



may be more complicated ( Mozurkewich (2000)| ). The "C2 and "'Cn schemes are called Pairwise 



and All-On-One combinations in Mozm'kewich (2000). 



Baseline redundancy is a factor that affects the SNR of a synthesis image for a given number 
of apertures. Since the number of large telescopes combined to form an interferometric array 
is likely to be limited, it is important to gain uv coverage by using nonredundant baselines. 
Throughout this paper, we only consider interferometers with nonredundant baselines. 

There are two types of noise which appear in aperture synthesis, one of which is removable 
while the other is intrinsic and not removable. Sampling noise is potentially introduced by the 
incompleteness of uv coverage. However, deconvolution techniques such as CLEAN and Maximum 
Entropy Method (MEM) ( Perley et al. (1989)| ), when applied to radio interferometric data, appear 



to compensate for sampling noise. We thus regard sampling noise as removable noise. On the 
other hand, noise resulting from the photodetection process is not removable. The SNR of a dirty 
image to which deconvolution has not been applied is limited by photodetection noise and detector 
read noise. Here we regard the SNR of a dirty image as the SNR of the synthesized image. 

There are two different methods by which a synthesized image can be constructed from 
the visibility data: inversion without total counts and true inversion. In the former, the zero 
frequency phasor is neglected, implying that the total number of photons in the synthesized 
image is zero. Despite this seemingly unattractive feature, this is the standard method in radio 
astronomy. In the latter method, the van Cittert-Zernike theorem is strictly applied and the zero 
frequency phasor as well as the n{n — l)/2 complex phasors is used. The image produced by this 
technique has the desirable property of nonnegativity. In the presence of additive background 
which overwhelms the signal, the magnitude of the zero frequency phasor is much greater than 
those of other phasors. The zero frequency phasor determines the DC offset level due to thermal 
background in the synthesized image. Not only is this background of little astronomical interest, 
it is also a major source of Poisson fluctuations. In this paper, we present the results obtained 
for both inversion with and without total counts. However, in practice faint source detection is 
limited by OH airglow background at J and H and thermal background in the thermal infrared. 
So we obtain those detection limits for inversion without total count. 

The paper is organized as follows. First we introduce two types of interferometers, the "(72 
interferometer and "C„ interferometer with expressions for SNRs of the images obtained by these 
interferometers in §2 and §3. We discuss limiting cases for spectrally broad light in §4. In §5, we 
present the sensitivities of interferometers composed of nine 10-m telescopes in comparison with a 
single 30-m telescope with adaptive optics. We discuss further on the sensitivity analysis in §6. 



2. "C2 interferometer 

Let there be n identical principal apertures from which we derive n main beams. Each beam 
is divided into n — 1 subbeams by the use of beam splitters. The resulting n{n — 1) subbeams are 



combined pairwise at rib = n[n — l)/2 detectors where rih is also the number of basehnes. The 
fringe pattern is formed at the focal plane of each detector which is either a one dimensional array 
of identical pixels or a two dimensional array of identical pixels which range from 1 to P (Figure 
1). For spectrally broad light, the fringe can be dispersed in the cross fringe direction with the 
use of a two dimensional detector. For a focal plane interferometer of this type, the influence of 



source shot noise on fringe phasor estimation has been fully formulated ( Walkup fc Goodman 



(1973) , Goodman (1985)] ). Let the interferometer be illuminated by a source and spatially smooth 



background. The intensity pattern on the rth detector (rth baseline) is given by 

< /r(x) >= 2 < lo > 1 + 7r COs(kX • B^/d + (t)r) , (l) 

where < /q > is the average intensity in each subbeam at the detector, B^ is the baseline 
vector, K = 27r/A is the wave number, and d is the distance between the aperture plane and the 
detector, x is the spatial vector in the detector plane, and 7j.exp(ii?i>r) is the complex visibility 
function at the baseline vector B^. In deriving (||), we have assumed that the incident light is 
spectrally narrow so that the fringe visibility depends on only the spatial correlation of the field. 

Because of the presence of background illumination, 

h = ^+ll (2) 

where /g and /q are source and background intensities respectively. Let 7^ be the fringe 
visibility of the source in the absence of the background, then 

-'0 "T-'o 

The complex visibility function 7^exp(i(/)r) is determined by the uv coordinates of the r th 
baseline (n, u) and the source brightness distribution on the sky S{x,y) by a Fourier transform 
relation: 

s rj,\ J S{x,y)exp{-2Tri{ux + vy)}dxdy 

j S{x,y)dxdy 

In an effort to reduce the clutter in the equations we henceforth drop the vector notation, but 
bear in mind that spatial frequencies, pixel locations, etc. are really vectors. The photoelectron 
detection theory ([Walkup fc Goodman (1973) ) takes into account the discrete nature of both 



photons and detector pixels. The average photoelectron count < kr{p) > at the pixel location 
specified by the integer index p of the detector is proportional to < I{x) >: 



< kr{p) >= 2 < Kq > 1 + 7r COs{pUJr + (f)r) ■ (5) 



Here, < ... > denotes averaging over the photoelectron-detection process. The product pujr is 
understood to be the scalar product of the pixel position vector p and the spatial frequency vector 
iOj. expressed in inverse pixel units. 

Let < C > be the average number of photoelectrons detected by the entire array in one 
integration period, and let 2< A^ > be the average number of photoelectrons per detector per 
integration time. Clearly then, < C >= 2 < N > Ub, and thus < N >=< C > /{n{n — 1)}. 
According to @, the average number of photoelectrons per detector is equal to 2 < Kq > P, and 
thus < Ko> P =< N >. 



2.1. Inversion without total count 

An analysis based on the photodetection theory gives the mean image /i as follows (NMOl). 



Ii{q) = < N >Y^ 7r cos{uJrq + (/)r) 

r 

-In 



r -'0 + -'0 



(6) 



where q refers to the pixels in the synthesized image; in particular q ranges from -Q/2 to 

+Q/2. 

The variance y[ii(g)] in the synthesized image ii{q) is given by 



T.r- / M r ^T Pa^, <C> Pa^ 

For a point source (7^ = 1) at the phase center [cpr = 0), the SNR of the synthesized image is 

/i(0) _(<C>/2){/oV(/o' + /o')} 



VVihm v/< C > /2 + UbPayi 



2.2. Inversion with total count 



(8) 



An analysis which takes into account shot noise and detector read noise gives the mean image 
I2 as follows 



hiq) = <N >J2 br cos {uJrq + <pr ) + 1] 



= <iV>^[— -^7rCosK(? + 0,) + l]. (9) 

r -'O + -'O 

The variance ^[^2(9)] in the synthesized image i2iQ) is given by 

V[l2{q)\ = rihi \ — )+ < N > 2_^'jrCOs{uJrq + (l)r) 

3<C> 3Pa^ .r V- 

= \-nh — \- < N > ^'yrCOS{UJrq + (Pr)- (10) 

r 

This variance is newly derived for this work. For a point source at the phase center, the SNR 
of the image is 

^2(0) _ <N>^,{2IS + l',)/m+l'o) ^^^^ 



\/yMO)\ y/(3 < iV > /2 + 3Pa2/4)+ < A^ > 1^/(1^ + 4) 



3. "C„ Interferometer 

In the "C„ interferometer, ah the n beams interfere on a single detector and n^ fringes are 
superposed. Both the baseline configuration and the detector can be either one dimensional or two 
dimensional. As an example, a two dimensional ^Cs interferometer is schematically presented in 
Figure 2. A special case is the combination of a one dimensional baseline configuration and a two 
dimensional detector for which the superposed fringes are dispersed in the cross fringe direction so 
that spectrally broad light can be used without bandwidth smearing (Figure 3). We call an "C„ 
interferometer of this type an "C^ interferometer. 

Let the interferometer be composed of n identical apertures and let it be illuminated by a 
source and spatially smooth background. The classical intensity distribution of the interference 
pattern by the n apertures has the average value 

< I(x) >=< /o > [n + 2j2 Igh cos(kx • Bgh/d + (t)gh)\ , (12) 

g<h 

where the various symbols have meanings similar to those in (||). gh denotes the baseline gh 
corresponding to the apertures g and h. Let < A;(p) > denote the photoelectron count distribution 
due to < /(x) >. As in §2, we discontinue the vector notation, assume that the total number of 
pixels is P, and note that < k{p) > is proportional to < I{x) >: 

< k{p) >=< Qo> n + ^-fgh cos{pujgh + (pgh) ■ (13) 

g<h 



Here < Qo > has approximately the same meaning as < Kq > in §2. However, since there is 
no beam spUtting, < Qq >= (n - 1) < Kq >■ We also define < M >= P < Qo >=< C > jn. 



3.1. Inversion without total count 

We find the mean synthesized image I3 to be 



h{q) = <M >J2lij cosiquij + (pij). (14) 

i<j 



The variance ^[^3(9)] in the synthesized image isiq) is given by 



n<M > +Pa^ _ , „, ^ /g _, 

rs _|_ jb 



VMq)] = 5 n6+<M>(n-2)^— -^7f,cos(g^,, + ,/.i,-). (15) 



This variance is derived in NMOl. For a point source at the phase center, the SNR of the 
image is 

/3(o) _ <M>{mm+i^)}^, 



V^[^3(0)J J{n<M> +Pct2)/2+ <M>{n- 2)m/m + I^)} 



<C> m/m + I^)}^in-l)/n 



<C> +Pa^+ <C> m/m + 4)}{2{n - 2)/n} 



(16) 



3.2. Inversion with total count 

We find the mean synthesized image I4 to be 



h{q) = < M > [| + ^ 7ij cos{quJij + (/.y)] (17) 

= — 5 \ 2^ 7ij cos{quJij + (pij). (18) 

Zi it.. 

The variance ^[^(9)] in the synthesized image iA{q) is given by 



V[H{q)] 



<M > 



n[nh 



+ 7^) + 2(n - 1) XI 2(J - l)7ii cos{quJij + 



i<j 



+Pa^ 



2 4 



(19) 
(20) 



This variance is newly derived for this work. For a point source at the phase center, the SNR 
of the image is 



/4(0) <C>n 

"TWM ~ x/<C > (3n^ - 5n + 3) + Pa'^{n^ - rTFI) 



(21) 



4. Limiting Cases for Spectrally Broad Light 

In this section, we consider the source shot noise hmit, background shot noise limit, and 
detector read noise limit for a point source and an extended source. For each case, the results of 
inversion with and without total count are shown. 

Formulae for means, variances, and SNRs were obtained for spectrally narrow light in the 
previous section. For a two-dimensional "Cm interferometer, the maximum spectral bandwidth is 
limited to about D/B where D is the aperture diameter and B is the maximum baseline length. If 
we use a bandwidth broader than this, fringes will be smeared out. In order to disperse fringes for 
a general two-dimensional baseline configuration, pupil remapping which transforms the exit pupil 
to a "C^-type interferometer is necessary. The effect of pupil remapping on the signal-to-noise 
ratio analysis is beyond the scope of this paper. However we can imagine that the SNR of the 
one-dimensional "C^ interferometer is a good measure of the SNR of an interferometer with pupil 
remapping. For "(72 and "C^ interferometers, fringes can be dispersed in the cross fringe direction. 
The formulae for SNRs obtained for spectrally narrow light hold also for dispersed fringes as can 
be easily verified as follows. 

We here consider the case with the "C2 interferometer and inversion without total count. Let 
Px be the number of pixels in the fringe direction and Py be that in the cross fringe direction. 
Therefore P = PxPy We assume that the photon spectrum within the bandwidth is not steep and 
the number of photons in one detector row is given by < C > /Py. Then the SNR of an image 
synthesized from a set of fringes in detector rows corresponding to ^/Py of the spectral bandwidth 
is given by 



snr 



_ < N > /Py J2r Ir COs{uJrq + (pr 

W(< N> /Py + P,a2/2) 



(22) 



If we further assume that the total bandwidth is narrow enough so that the source structure 
does not change within it. In this case, Py images can be coadded to form a final image. The SNR 
of the final image is given by 



SNR = ^|PySnr 



< N >J2r^r COs{cL!rq + 4>r) 

(23) 



/n6(< N > +P^Pya^/2) 

< N > J2r 7r COs{LJrq + <Ar 



^nb{< N > +Pay2) 

This happens because the photon count < A^ > or < C > and the pixel number P are 
linear in the expressions of the mean and variance. The situation is the same for the other SNR 
expressions. 



4.1. Numbers of Pixels 

Let the central frequency be u and the bandwidth be Az/. Then there are 2iy/Aiy fringe cycles 
in each interferogram. For the Nyquist sampling theorem, Ai^/Ai^ pixels are needed in the fringe 
direction. 

For the "C2 and "C^ interferometers, we consider the number of pixels in the cross fringe 
direction (dispersion direction) in each dispersed fringe. The fractional bandwidth Av* jv for which 
bandwidth smearing is small is given by Av* jv < D/B. If we require that Av* jv = D/{2B), 
bandwidth smearing will be negligible. For a spectral bandwidth of Az/, Ai'/A*iy = 2{Aiy)/(D/B) 
pixels are needed as the number of pixels in the cross fringe direction. So the number of pixels of 
the two dimensional detector is P = Au/Au x 2{Av)/ {D/B) = 8{B/D). 

For the "C„ interferometer, the longest baseline can be aligned to one of the two sides of a 
two dimensional detector and the minimum number of pixels is P = {Au/Au)'^ = 16{B/D)'^ and 
the maximum allowed fractional bandwidth is D/B. 



4.2. SNR Tables 

SNR expressions of "C2, "C„, and "C^ interferometers are given for a point source and an 
extended source, for source shot noise limit, background shot noise limit, and detector read noise 
limit, and for inversion with and without total count in Tables 1 through 5. In some cases the 
SNR expressions are not well defined and omitted from the tables. 

We define a quantity F^ which is proportional to the photon spectrum of the source, collecting 



10 



area, and integration time. We express the total photon count for the "C2 and "C^ interferometers 
as < C >= Fy/S.u. Then the total photon count for the two-dimensional "C„ interferometer is 
< C >= Fyv{D /B). The SNRs are expressed in terms of Fy, Az^, and D /B. 



4.3. Point Source/Source Shot Noise Limit 

The case with a point source in the shot noise limit is given in Table 1. 
By taking ratios of the SNR expressions, we find 



^-./^™- /^ J|^. (-) 



/ 3„ _ 4 
SNR,/SNE!, = J,^—^>1, (25) 



s.n,s.u..^-^^l^^^-^. <., 



SNR.ISNR', = ./ '"='"' :f + '' > 1, (27) 

V 5n^ 



SNR2/SNRi = sl-. (28) 

In a typical situation, ^v jv > D jB and the "C2 and '^C^ interferometers are superior to 
the two-dimensional "C„ interferometer. As the expression (^) shows, the "(72 interferometer is 
superior to the "C2 interferometer by a factor ^12/5 in SNR for a large n for inversion with total 
count. In the source shot noise limit, the total count information is genuinely of the source. The 
highest SNR is obtained for the ""(72 interferometer for inversion with total count (eq. (E^ 



4.4. Background Limit 

In the background limit, only inversion without total count is the valid method of analysis, 
since the total count is dominated by background photon counts. The SNR expressions are given 
in Table 2. 

By taking the ratios, we find 
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SNR,/SNR, = r^J^, (29) 



and 



SNR,/SNR'^ = J^^^ < 1- (30) 

As in the case with the source shot noise limit, the two-dimensional "C„ interferometer has a 
disadvantage of the narrow bandwidth. For a large n, the "C^ interferometer is slightly superior 
to the "C2 interferometer by a factor \/2 in SNR. 



4.5. Read noise limit 

In the read noise limit, the numbers of pixels for the "C2 and "(7^ interferometers are 
P = 8{B/D) and that for the two-dimensional ^Cn interferometer is P = 16{B / D)^ . The SNR 
expressions are given in Table 3. 

If we take the ratios of the SNRs, we find 

SNR,/SNRs = -^ — {^f'\ (31) 

n — \ u D 

SNRi/SNR'^ = ^— < 1, (32) 



SNR.ISNR, = 1 &^^(§)=/^ (33) 

n V 3 n^ — n) V D 



SNR^ISNg, = l{^^ < 1. (34) 

The ^C'^ interferometer is always superior to the "C2 interferometer because of the smaller 
number of detector pixels. The two-dimensional ""Cn interferometer has a limited use for a small 
B/D. 
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4.6. Extended Source/Source Shot Noise Limit 

For a fully extended source, there is no signal other than the total count. Only inversion with 
total count is the valid method of analysis for the fully extended source. The SNR expressions are 
given in Table 4. 

By taking the ratios, we find 



SNR,/SNR, = f^^^^^^. (35) 



SNR2/SNR'^ = J!!l_^1±1 > 1. (36) 

The "C2 interferometer is superior to the "C„ and "C^ interferometers. 



4.7. Extended Source/Background Limit 

For an extended source, there is no signal other than the source total count. However, the 
total count is dominated by background photons. Therefore neither inversion without total count 
nor inversion with total count is the valid method of analysis. 



4.8. Extended Source/Read Noise Limit 

As in the case with extended source/shot noise limit, the only valid method of analysis is 
inversion with total count. The SNR expressions are given in Table 5. 

By taking ratios, we find 



V 3 n^ — n) V D 



and 



SNR^ISNR!^ = / 2("^ /^+1) < 1. (38) 

y 3(n^ — n) 

The two-dimensional "C„ interferometer is useful only for a small B /D. For large n, 
SNR2/ SNR!^ attains an asymptotic value of a/2/3. 
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5. Sensitivities of Interferometers Composed of Nine 10-m telescopes 

In this section, we calculate the sensitivities of "C2 and "C^ interferometers composed of nine 
10-m telescopes and compare them with the sensitivity of a 30-m single telescope with adaptive 
optics. lO-o" detection limits for a point source are obtained for one-hour integration and plotted 
in Figure 4. 



5.1. J and H 

Throughput of the interferometers is assumed to be 10%, while that of the single telescope 
with adaptive optics is assumed to be 60%. The read noise level is assumed to be 10 electrons 
per pixel and the maximum baseline length which determines the number of detector pixels is 
assumed to be 200 m. The coherent integration time is assumed to be one second, while a larger 
value will decrease the effect of detector read noise. 

It turned out that both read noise and sky background due to OH airglow affect the 
performance of the interferometers. We evaluate the SNRs for inversion without total count 
because sky background is much greater than the signal near the detection limit. The sensitivity 
of the telescope is limited by sky background. At J, limiting magnitudes are 24.7, 25.7, and 28.0 
respectivity for the "(72 interferometer, "C^ interferometer, and the single telescope. At H, they 
are 23.6, 24.0, and 26.2 respectively. 



5.2. Thermal Infrared Wavelengths 

For the interferometers, we evaluate the SNRs for inversion without total count because 
the instrument background is much stronger than the source signal near the detection limits. 
The temperature of interferometer optics is assumed to be 273 K and the optical throughput is 
assumed to be 10%. Therefore the background emission is estimated from a 273K blackbody with 
90% emissivity. For the single telescope with adaptive optics, the throughput of the optics is 
assumed to be 60%. The temperature of the optics is assumed to be 273 K, and its emissivity is 
assumed to be 40%. 

Longward of 2 ^um, the sensitivities do not depend on beam combination geometry very 
much. This implies that detector read noise is negligible. At 2.2 /xm, the sensitivities of the 
interferometers are around 330 nJy or 22.9 mag, while that of the single telescope is 10 nJy or 26.9 
mag. Longward of 5 //m, the single telescope is 5.8 mag more sensitive than the interferometers. 
At 10.5 fim, the sensitivities of the interferometers are about 2 mJy or 10.7 mag, while that 
of the single telescope is 57 /iJy or 14.4 mag. At wavelengths shorter than 2 ^um, the baseline 
length affects the sensitivities of the interferometers due to the change in the number of pixels, 
while longward of 2 /im, it does not matter as long as the source is unresolved. Here we do not 
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consider the case for a fully extended objects. As we have seen from the previous section, the 
interferometers are not sensitive to fully extended objects under background limited conditions. 



5.3. Loss-Free Interferometers 

In the above evaluation of sensitivities, it is assumed that the throughput of the interferometers 
is 10%, while that of the single telescope is 60%. It is natural to question whether the difference 
in sensitivities is due to that in throughput or something more fundamental. 

As discussed by PK89, there is a basic difference between an interferometer and a single 
telescope in that the former takes the correlation of signal from multiple apertures. However the 
loss of SNR by this process is only a factor of v2 in the case of an ideal "C2 interferometer with 
inversion without total count in the source shot noise limit. A similar situation is expected in the 
case of background shot noise limit where the number of detector pixels does not affect the SNR. 
Then the only situation in which the difference is conspicuous is the read noise limit where the 
interferometer requires many more detector pixels to sample fringes. It is of theoretical interest to 
assume ideal interferometers and compare their performance with that of the single telescope. 

Again we compare the interferometers with nine 10-m telescopes and a 30-m single telescope. 
Since each of the 10-m telescope and the 30-m single telescope require adaptive optics, the 
throughputs are assumed to be 60%. The sensitivity limits are plotted in Figure 5, where the 
difference is about by a factor of five. Therefore the fundamental difference between in the 
interferometers and the single telescope is rather small and the large difference seen in the previous 
analysis is due to that in throughput. 



5.4. Baseline Dependence 

The SNR expressions of the interferometers are not explicitly baseline dependent. Since 
the astronomical source is assumed to be unresolved, only the total number of detector pixels is 
dependent on the lengths of baselines. Detector read noise is more significant as the baselines 
become longer. In order to use long baselines, it is essential to keep detector read noise low so 
that the SNR is not limited by read noise. 

For a source of finite extent, long baselines start to resolve it and the visibilities 7^ of those 
baselines decrease from unity and the phases (pr become less coherent depending on the symmetry 
of the source. For the extended source, the SNRs of the interferometers are calculated as follows. 
For a particular baseline configuration, the complex visibilities are obtained by the Fourier 
transform relation (4) and the number of detector pixels are calculated as above. By inserting 
those numbers and taking the summation over different baselines, one can obtain the SNRs of the 
dirty images. 
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6. Discussion 

6.1. Interferometry for extragalactic astronomy 

One of the motivations of tliis analysis is to evaluate the capability of interferometry for 
extragalactic astronomy. We inevitably have to consider extended objects for this purpose and take 
into account the guide star availability. Fringe tracking requires the presence of a natural guide 
star in the infrared. If the adaptive optics also needs the same natural guide star, the limiting 
magitude is around R = 15, and an isoplanatic angle is order of 10" at 2 /xm. Requirements for a 
natural guide star limit the range of observing targets to bright active galactic nuclei (AGNs) and 
bright quasars. 

The central continuum source of an AGN or a quasar will act as a natural guide star and the 
fringe visibilities are close to unity and phases are close to zero. Therefore SNR calculations made 
in the previous section are applicable to this case. Small deviations in a reconstructed image from 
a point source indicate the presence of a faint extended structure such as jet. The 10a limiting 
magnitude of K=25 indicates SNR = 10^ for a K=15 quasar. In other words, jet can be searched 
for with a dynamic range of 10^:1. 

Now we consider resolving broad-line regions (BLRs). According to the photoionization model 
( Davidson fc Netzer (1979) ), the size of a BLR scales as \/L where L is the luminosity of the AGN. 



For an AGN or a nearby quasar toward which the geometry is Euclidean, the observed flux / is 
given by / = L/d? where d is the distance to the object. Then the angular size oc ^/f. Therefore 
the brightest quasar has the largest BLR in angular size. 3C273 is the brightest quasar (K=10) 
and easiest to observe because Paa is redshifted to K band. The necessary spatial resolution is 
order of 1 milliarcsecond or the necessary baseline length is 400 m. The FWHM of the line is 
3400 kms~^ and the line/continuum ratio below 8000 kms~^ is 0.08 ( Sellgren et al. (1983)| ). In 



order to obtain the spatial structure of the BLR, we wish to spectrally resolve the broad line or at 
least separate the blue and red components. If we require a spectral resolution of 3000 kms"^, the 
continuum flux per spectral resolution is K=14 equivalent. A spectral image will be composed of 
an unresolved source of K=14 and a resolved source with 8% of the flux of the unresolved source. 
Again, the complex visibilities are close to those of a point source and the SNR argument made 
above for jet applies to this case as well. We expect that the BLR of 3C273 will be resolved with 
a dynamic range more than 10^:1. At least for AGNs and quasars, infrared interferometry with 
large telescopes appears effective. 



6.2. Comments on RR99 

In this subsection, we compare our results with that of RR99 on Fizeau-type interferometers 
or two-dimensional ""C^ interferometers. First of all, for spectrally broad light, we have seen that 
two-dimensional "C„ interferometers or Fizeau-type interferometers have limitations in the ratio 



16 



of the telescope aperture diameter and the basehne length to avoid bandwidth smearing. This 
is not discussed in RR99. Second, RR99 states that the "C„ interferometer is superior to the 
"(72 interferometer. However, we have seen from the analysis of limiting cases that the situation 
is not so simple. As PK89 noted, there is no major difference in performance between both 
interferometers in the source shot noise limit. This is also true with the case for the background 
limit. In the read noise limit, the "C^ interferometer is superior to the "C2 interferometer due to 
the difference in the number of detector pixels. 

These differences originates from that in the definitions of the SNRs. RR99 defines the SNR 
as that of the central peak in the interferogram directly formed by beam combination. On the 
other hand, we define the SNR as that of the dirty image formed by aperture synthesis. In the 
latter, the interferogram formed by beam combination is Fourier transformed to measure complex 
visibilities which fill the uv plane and then they are Fourier transformed back to the image plane 
to form a dirty image. In RR99, the detector pixels of concern at the focal plane are only a 
few pixels at which the peak is located, while in the latter many pixels which sample the entire 
interferogram are considered. 



6.3. Remaining Issues 

We have obtained the sensitivity limits of the interferometers for perfect wavefront correction 
by adaptive optics and perfect fringe tracking. In reality, both adaptive optics and fringe tracking 
require a guide star and there is a limitation in the precision of wavefront and delay compensations. 
Imperfection in phase directly influences the fringe visibility 7 which is proportional to the SNRs. 
An analysis which takes into account this imperfection in phase is difficult because introduction of 
atmospheric disturbances opens up a vast parameter space. Our analysis corresponds only to the 
best case designing. 

In this paper, we considered only nonredundant baselines. However, if we consider a closely 
packed array of telescopes, we have to take into account the redundancy of the baselines. Although 
the behavior of SNRs are expected to be the same as the nonredundant case, deriving formulae of 
SNRs for the redundant case is non trivial or analytically impossible. 

Finally one essential and difficult question is whether the definition of the SNR as that of the 
dirty image is valid. In other words, we have left out the issue of deconvolution by assuming that 
it is perfect. We believe that this is practically of no problem because the SNR of the dirty beam 
is usually much higher than that of the dirty image. Nevertheless deconvolution is a nonlinear 
process and the error propagation through it is not fully understood. 



17 



7. Conclusions 

We present the formulae of SNRs of interferometers for the purpose of comparing ground-based 
interferometers composed of multiple large telescopes with a single giant telescope with adaptive 
optics. In deriving sensitivity limits, perfect wavefront correction by adaptive optics and perfect 
cophasing by fringe tracking are assumed. We consider the case in which n beams are pairwise 
combined at n{n — l)/2 detectors ("C2 interferometer) and the case in which all the n beams are 
combined at a single detector ("C^ interferometer). 

(1) In the read noise limit, the "Cn interferometer with one dimensional baseline configuration 
is superior to the "C2 interferometer, while in the background limit, the advantage of the "C„ 
interferometer is small. 

(2) We compare the point-source sensitivities of interferometers composed of nine 10-m diameter 
telescopes and a 30-m diameter single telescope with adaptive optics between 1 and 10 ^m for 
lO-cr detection in one hour. Between 1 and 2 /xm, the sensitivities of the interferometers are 
limited by detector read noise and sky background, while longward of 2 /im, they are limited by 
the instrumental thermal background. At K band, the sensitivities of the interferometers are 
around 23 mag, while that of the single telescope is 27 mag. 

I thank H. Matsuhara for allowing me to use figures and the referee P. Lawson for useful 
comments. 
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Fig. 1. — Schematic diagram of a focal plane interferometer. There is one detector for each baseline. 

Fig. 2. — Schematic diagram of a ^Cs interferometer. Beams from all apertures are superposed on 
a single detector. 

Fig. 3. — Schematic diagram of a one-dimensional ^Ca interferometer. Fringes are spectrally 
dispersed in cross fringe direction. 

Fig. 4. — Detection limits of the "C2- and ID "C^ interferometers and that of the single giant 
telescope are plotted as functions of wavelength. Symbols are sensitivities are at J and H of "C2 
(asterisk), ID "C„ (square) and telescope (circle). Assumed parameters are described in the text. 

Fig. 5. — Loss-Free Interferometers. Throughputs of the interferometers and the single telescope 
are assumed to be the same. Meanings of the symbols are the same as in Figure 4. 
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Table 1. SNR (Point Source / Shot Noise Limit) 



Case 



Notation 



Expression 



"C2/inversion w/out total count SNRi 

"C2 /inversion with total count SNR2 

"C„ /inversion w/out total count SNRs 

"Cn/ inversion with total count SNR4 

"C^/inversion w/out total count SNR^j 

"C^/ inversion with total count SNR'4^ 



V^: 



TTTKv 



IF^^Ku 



^^\/FvV% 



Zn^ — 5n+3 



/fu 



i^/KKt 



3n^ — 5n+3 



VKAv 



Table 2. SNR (Point Source / Background Limit) 



Case 



Notation 



Expression 



"C2/ inversion w/out total count SNRi 

"Cn/ inversion w/out total count SNR3 



^C^/ inversion w/out total count 



SNm 



'0 

pn 






-i^i^Uf 



/^/fUa^ 



Table 3. SNR (Point Source / Read Noise Limit) 



Case 



Notation 



Expression 



"C2 /inversion w/out total count SNRi 

"C2 /inversion with total count SNR2 

"C„/inversion w/out total count SNR3 

"Cn/inversion with total count SNR4 

"C^ /inversion w/out total count SNR-j 

"C^/ inversion with total count SNR'^ 



4,y/nt(B/D) <^ 

1 fVAi/ 

76 nt,(B/ D) <^ 

4 y n (b ) CT 
i / n^ 'tD\2 F„v 
4 y n'-'-n+l ^ B' ct 
1 /n-1 Fi,Av 



1 

VkbJd) 



n^ —n-\-l <y 
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Table 4. SNR (Extended Source/Source Shot Noise Limit) 



Case 



Notation Expression 



"C2/inversion with total count SNR2 

"C„/inversion with total count SNR4 
"C4/inversion with total count SNR'4 



V^ 



F,yAU 



F^v{D/B) 



-n+1 



F,^Ai/ 



-n+1 



Table 5. SNR (Extended Source/Read Noise Limit) 



Case 



SNR 



Expression 



^6*2 /inversion with total count SNR2 
^Cn /inversion with total count SNR4 
^C^/inversion with total count SNR4 







1 


F,.l 


\v 




yi2l 


[n^-n)(B/D) 
1 


(7 


(B/D) 


Vl6(n2- 


-n+l){S/D)2 
1 


F„ 


a 



^8{n2-n+l){S/D) 



This figure "figl.gif ' is available in "gif ' format from: 



http://arXiv.org/ps/astro-ph/0107262vl 



This figure "fig2.gif" is available in "gif ' format from: 



http://arXiv.org/ps/astro-ph/0107262vl 



This figure "fig3.gif" is available in "gif ' format from: 



http://arXiv.org/ps/astro-ph/0107262vl 
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